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To avoid the complicated topology of surviving clusters induced by standard Strong Disorder 
RG in dimension d > 1, we introduce a modified procedure called 'Boundary Strong Disorder RG' 
where the order of decimations is chosen a priori. We apply numerically this modified procedure 
to the Random Transverse Field Ising model in dimension d — 2. We find that the location of the 
critical point, the activated exponent if) ~ 0.5 of the Infinite Disorder scaling, and the finite-size 
correlation exponent vp$ — 1.3 are compatible with the values obtained previously by standard 
Strong Disorder RGOur conclusion is thus that Strong Disorder RG is very robust with respect to 
changes in the order of decimations. In addition, we analyze in more details the RG flows within 
the two phases to show explicitly the presence of various correlation length exponents : we measure 
the typical correlation exponent v tyv — 0.64 in the disordered phase (this value is very close to the 
correlation exponent Vp Ure (d = 2) ~ 0.63 of the pure two-dimensional quantum Ising Model), and 
the typical exponent Uh — 1 within the ordered phase. These values satisfy the relations between 
critical exponents imposed by the expected finite-size scaling properties at Infinite Disorder critical 
points. Within the disordered phase, we also measure the fluctuation exponent oj ~ 0.35 which is 
compatible with the Directed Polymer exponent cjdp(1 + 1) = | in (1 + 1) dimensions. 



I. INTRODUCTION 

Strong disorder Renormalization (see [l| for a review) has been first introduced for one-dimensional quantum spin 
chains [2H4J , where exact solutions can be obtained because the renormalized lattice of surviving degrees of freedom re- 
mains one-dimensional. In dimension d > 1, the Strong Disorder RG procedure cannot be solved analytically, because 
the topology of the lattice changes upon renormalization, but it has been implemented numerically, in particular for 
the quantum Ising model [5l-[l5l]. Nevertheless, the complicated topology that emerges between renormalized degrees 
of freedom in dimension d > 1 tends to obscure the physics and slow down the numerics, because a large number of 
very weak bonds are generated during the RG, that will eventually not be important for the forthcoming RG steps. 
Various types of simplifications have been thus proposed, like the 'maximum rule' [5rJi"oj possibly supplemented by 
some very efficient algorithm p^4l5j | , the introduction of a cut-off within the full sum rules [l6| or the planar approx- 
imation [17j | . Recently we have proposed to follow another strategy : the idea is to allow some changes in the order of 
decimations with respect to the full procedure in order to maintain a simple spatial renormalized structure. We have 
already applied this idea in two ways : (i) in fl8| . we have proposed to include strong disorder RG ideas within the 
more traditional fixed-length-scale real space RG framework that preserves the topology upon renormalization, with 
numerical results for various types of fractal lattices; (ii) in jl9[, we have proposed for the Cayley tree geometry some 
'Boundary Strong Disorder RG procedure' that preserves the tree structure, so that one can write simple recursions 
with respect to the number of generations. In both cases, we have checked that in dimension d = 1, these modified 
procedures correctly capture all critical exponents except for the magnetic exponent [3 which is related to persistence 
properties of the full RG flow. In the present paper, we adapt this idea of 'Boundary Strong Disorder RG procedure' 
to the two-dimensional case and present the corresponding numerical results, that we compare with the results of 
standard Strong Disorder RG [5l4l5| and with quantum Monte-Carlo [2fl [2l|. 

The paper is organized as follows. In section UH we define the Boundary Strong Disorder RG procedure for the 
two-dimensional square lattice. In the following sections, we discuss the numerical results obtained by this procedure. 
In the disordered phase (section ITTT|) . we measure the typical correlation exponent v tyv , the fluctuation exponent to and 
the essential singularity exponent k. In the ordered phase (section IIV|) . we measure the typical correlation exponent 
Vh- In the critical region (section [Vj , we find that the location of the critical point, the activated exponent tp and the 
finite-size correlation exponent vps are compatible with the values obtained previously by standard Strong Disorder 
RG [5l-[l5j. Our conclusions are summarized in section ["VTl 

II. BOUNDARY STRONG DISORDER RG PROCEDURE IN d = 2 

As recalled in Appendix the Strong Disorder Renormalization for the quantum Ising model is an energy-based 
RG, where the strongest ferromagnetic bond or the strongest transverse field is iteratively eliminated. In this section, 
we introduced a modified procedure, called Boundary Strong Disorder RG, that preserves a simple spatial structure. 
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A. 



Initial model 



In this paper, wc consider the quantum Ising model defined in terms of Pauli matrices 



n = - 




.X 



(1) 



<i,j> 



on the square lattice in dimension d — 2 where the initial nearest-neighbor couplings J™ are independent random 
variables drawn with the box distribution on the unit interval [0, 1] 



nj(J^) = 6(0 < < 1) 



(2) 



and where the initial transverse fields h] m > are independent random variables drawn with the box distribution on 
the interval [0, h] 



so that the parameter h is the control parameter of the quantum phase transition as in Refs [7H15L I20L l21j . 

For the numerical results, we consider more precisely a square lattice containing (2L — l) 2 spins of coordinates 
(x = 1, 2, .., 2L — l;y = 1, 2, .., 2L — 1. Each spin has its random initial transverse field drawn with the distribution 
of Eq. [3] and is connected to four neighbors via random ferromagnetic couplings drawn with the distribution of Eq. 
[3J All exterior sites situated along the boundaries at x = 0, x — 2L or y = 0, y — 2L are identified to a single formal 
'external spin' to keep track of the coupling to the boundary of the finite sample. 



We wish to eliminate sites in a simple deterministic order, starting from the boundary : we will first eliminate sites 
that are at distance 1 from the boundary having coordinates (x — 1) or (x = 2L — 1) or (y = 1) or (y = 2L — 1); then 
sites that are at distance 2 from the boundary having coordinates (x — 2) or (x = 2L — 2) or (y = 2) or (y = 2L — 2); 
and so on, up to sites that are at distance (L — 1) from the boundary having (x = L — l) or (x = L + 1) or (y = L — l) or 
(y = L + l). At a given stage of the RG, we have a renormalized spatial structure containing a 'corona of renormalized 
boundary sites' (the number of sites in the corona scales as the surface (L — l)^ 1 when the corona is at distance I 
from the boundary) and the 'interior of the corona' (the number of sites in the corona scales as the volume (L — l) d 
when the corona is at distance I from the boundary). We have the following properties : 

- the 'interior of the corona' contains sites that have not yet been modified with respect to the initial model, i.e. 
the sites are characterized by their initial random fields h lm (i), and are connected to their initial neighbors by their 
initial ferromagnetic coupling J™. 

- the 'corona' contains renormalized boundary sites (i) that have renormalized transverse fields h(i) and that are 
connected to the formal external spin via some renormalized coupling J ext (i). These corona sites are connected 
to 'interior spins' via their initial ferromagnetic coupling J*™. Finally, there may exist renormalized ferromagnetic 
couplings Jij between any two pair of sites belonging the corona. 

- the sites outside the 'corona' have been already eliminated, and these eliminations are responsible of the renor- 
malized variables characterizing the corona sites. 



To eliminate a given site i of the corona, we determine the maximum between its renormalized transverse field h(i) 
and its renormalized ferromagnetic couplings Jy with the other sites of the corona or interior sites 



MM?') = -^(0 < K' 1 < h) 



(3) 



B. Boundary Strong Disorder RG spatial structure 



C. 



Boundary Strong Disorder RG rules for a corona site 



Hi = max [hi, Jy] 



(4) 



(Note that the renormalized external coupling J ext (i) is excluded, since it is just a 'passive' variable used to measure 
the effective coupling to the initial boundary) 

Then we apply the Strong Disorder RG rules as follows (see Appendix |Aj : 
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i) If fli = hi, then the site i is decimated, and all couples (j, k) of neighbors of i are now linked via the renormalized 
ferromagnetic coupling 

jnew _ j , JjiJik ( r\ 
J 3 k — J ]k H [?) 

Accordingly, the external couplings of all neighbors j of i are renormalized according to 

j jext 

je Xt ,ne W = je X t + (g) 

J J K 

ii) If Vti — Jij, then the site i is merged with the site j. The new renormalized site j has a reduced renormalized 
transverse field 

h new = Ml ( 7 ) 
Jij 

This renormalized site j is connected to other sites k via the renormalized couplings 

J™ k ew = J jk + Jik (8) 
In particular, the external coupling of the renormalized site j becomes 

= jf * + jf xt (9) 

Note that when the site i is eliminated, all interior sites that were connected to the site i become sites of the new 
corona. 

In the final state of the RG procedure, only the center site of coordinates (x c = L,y c = L) remains. The observables 
under interest are its final renormalized transverse-field h^ and its ferromagnetic coupling J£ xt to the formal 'external 
spin', i.e. to the initial boundary of the square sample. To simply the notations, Jf^ 1 will be denoted simply by Jl 
from now on. In the following, we will concentrate on the typical values h 1 ^ 9 and </^ p defined by 

In h l ^ p = In hL 

\nfy p = ToJl (10) 



and on the widths of the distribution of In h^ and In Jl defined by 

1/2 



Mn/ir 



e ((ln/i L )2-(ln/i L ) 2 ) 
Am,;, ee ((lnJ i )2-(hTjL) 2 ) 1/2 (11) 



where the overbar denotes an average over the disordered samples. 



D. Numerical details 



We have followed numerically the Boundary RG rules for square samples containing Nl = (2L — l) 2 spins, for 
various sizes L < 100 corresponding to Nl < 9801. For a given size L, the number n s (L) of independent disordered 
samples we have been able to study depends on the value h of the initial disorder distribution of Eq. [3] So let us 
give some typical values we have used in the critical region (we were able to study more samples in the ordered phase 
h < h c and less samples in the disordered phase) 

L = 6,10,20,40,60,80,100 
n s (L) = 10 8 ,2.10 7 ,3.10 5 ,2.10 4 ,2.10 3 ,5.10 2 ,2.10 2 (12) 

Our various data shown below are compatible with a critical point located around the value (see the definition of the 
control parameter h in the initial disorder distribution of Eq. [3]) 

h c ~ 5.15 (13) 

which is sligthly lower but close to the values found previously u sing the standard Strong Disorder RG rules with 
the maximum rule, namely h c ~ 5.3 0] and h c ~ 5.35 [Tfj| - [T2i [Til ITsT ]. Let us first discuss the properties of the two 
phases, before we turn to the critical region. 
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III. ANALYSIS OF THE DISORDERED PHASE 



A. RG flow of the renormalized external coupling Jl in the disordered phase 



typ 

ln(-ln J ) 




In L 




FIG. 1: (a) Statistics of the logarithm of the external renormalized coupling (In Jl) in the disordered phase (here h = 8) : the 
RG flows of the typical value In j£ yp and of the width Ai n j L in a log-log plot display slopes 1 and u ~ 0.35 respectively (Eq. 
I19[) (b) Inverse of the typical correlation length £t yP of Eq. [15] as a function of the control parameter h > h c : the correlation 
length exponent v tyv of Eq. \W\ is of order v typ — 0.64 (see the log- log plot in inset) . ] 



In the disordered phase, the renormalized external coupling Jl is expected to present the following scaling 

lnJ L = — — + L u A{h)u 
styp 



(14) 



The first non-random term describing the exponential decay with the size L ( see Fig. [T] (a) ) defines the typical 
correlation length ^typ 



In J L yp = In Jl oc — - — 



(15) 



On Fig. [T] (b), we show how 1/^typ varies as a function of the control parameter h of the transition : our data are 
compatible with the power-law divergence (see the log-log plot in the inset of Fig. [1] (b)) 



ityp oc (A - h c y 



with a typical correlation exponent of order 



Vtyp 



0.64 



(16) 



(17) 



To the best of our knowledge, this is the first numerical measure of this typical exponent v ty p within the disordered 
phase, since previous studies have concentrated on the critical region where finite-size effects are governed by another 
correlation length exponent vps (see section [Vj) . We note that the value found here for the typical exponent of Eq. 
[T7] turns out to be very close to the correlation exponent Vp Ure {d = 2) ~ 0.63 of the pure two-dimensional quantum 
Ising Model . (The latter is known to coincide with the correlation exponent Vp^(d + 1 = 3) ~ 0.63 of the pure 
three-dimensional classical Ising Model as a consequence of the quantum-classical correspondence (23[). Since in 
dimension d = 1, the typical exponent v typ {d = 1) = 1 also coincides with the correlation exponent t , ^ ire (rf = 1) = 1 
of the pure one- dimensional quantum Ising Model (and equivalently with the exponent v^*£(d + 1 = 2) = 1 of the 
pure two-dimensional classical Ising Model ), it would be interesting to determine whether these coincidences continue 
in higher dimensions d > 3, i.e. whether the typical exponent takes the simple value v ty p{d > 3) = Vp Ure (d > 3) = 



<°/ e s (d+l>4) = | ? 

The second term in Eq [TJ] contains an 0(1) random variable u, which is expected to be subleading with respect to 
the first term, i.e. the width Ai„ j L of the distribution of In Jl is of order L u with some fluctuation exponent u) < 1 



Mn J L 



((lnJ L ) 2 -(lnJ L )' 



1/2 



L— >+oo 



(18) 
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We have argued in [22[ that this exponent ui should coincide with the droplet exponent uppiD = d— 1) of the Directed 
Polymer with D = {d — 1) transverse directions. For our present case in d = 2, the droplet exponent of the Directed 
Polymer is exactly known to be ujdp{D = 1) = | 

u = ujdp(D = 1) = - (19) 

in agreement with our numerical results shown on Fig. [T] (a). Again, to the best of our knowledge, this fluctuation 
exponent uj within the disordered phase had not been measured yet, since previous studies have concentrated on the 
critical region. 



B. RG flow of the renormalized transverse field hi, in the disordered phase 




FIG. 2: Statistics of the asymptotic finite renormalized transverse field /i.l->oo in the disordered phase as a function of the 
control parameter h > h c : the critical exponent k of Eq. [22] is of order k, ~ 0.65 (see the log-log plot in inset). 



In the disordered phase, the renormalized transverse held /il is expected to remain a finite random variable as 
L — > +oo. In particular, the typical value remains finite 



In h l £ p ee In h L oc In ^ (20) 

To analyze more clearly the statistics of the finite renormalization with respect to the initial random fields hi n i drawn 
with the box distribution of Eq. [3] corresponding to 

h dh- ■ 



Inhini = / — r — hi/iwM = hi/i - 1 (21) 



we show on Fig. [2]the difference (hihoc — hihi n i) as a function of h > h c . Our data are compatible with the expected 
essential singularity 



In^L oc -(h-h c )~ K (22) 

with an exponent of order 

k ~ 0.65 (23) 
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FIG. 3: (a) Statistics of the logarithm of the renormalized transverse field (In /il) in the ordered phase (here h — 1) : the RG 
flows of the typical value (ln/i^ /p ) and of the width A\ n h L in a log-log plot display slope d — 2 and slope d/2 = 1 respectively, 
(b) Inverse of the correlation length of Eq. [21] as a function of the control parameter h < h c : the correlation length exponent 
vh of Eq. [16] is of order Vh — 1. 



IV. ANALYSIS OF THE ORDERED PHASE 



A. RG flow of the renormalized transverse field III in the ordered phase 



In the ordered phase, the logarithm of the renormalized transverse field is expected to behave extensively in the 
volume L d 



ln^^ln/ii r oc -(^) (24) 



L 

with d = 2, in agreement with our data shown on Fig. |3J(a). The length scale £/, represents the characteristic size of 
finite disordered clusters within this ordered phase. It is expected to diverge as a power-law near the transition 

& oc (h c - h)~" h (25) 

The corresponding correlation exponent Vh plays in the ordered phase a role similar to v tyv in the disordered phase 
(Eq. IT5|) . Our data are compatible with a value of order (see Fig. [3] (b)) 

v h ~ 1 (26) 

Again, to the best of our knowledge, this is the first numerical measure of the typical exponent Vh within the ordered 
phase, since previous studies have concentrated on the critical region where finite-size effects are governed by another 
correlation length exponent vps ( see section [VJ. We note that this exponent takes also the same value Vh(d = 1) = 1 
in d = 1, but we are not aware of any argument in favor of this simple constant value as d varies. 

As shown on Fig. [3] (a), the width of the distribution of the logarithm of the renormalized transverse field grows 
linearly in L 

A lnhL = ((lnh L ) 2 - (hJ^) 2 ) 1/2 ex L (27) 



B. RG flow of the renormalized external coupling Jl in the ordered phase 



In dimension d — 1 where there is no underlying classical ferromagnetic transition, the typical renormalized coupling 
remains finite in the ordered phase, and presents the same essential singularity as in Eq. 1221 However in dimension 
d > 1 where there exists an underlying classical ferromagnetic transition, the renormalized couplings Jl is expected 
to grow at large L with the scaling of the classical random ferromagnetic model (see Fig. @] ) 

\nJ L vp =EM oc ]n(a(h)L d °) (28) 

L— >-\-oo 




- 



In A, 



■s— h-b- 



FIG. 4: Statistics of the renormalized external coupling Jl in the ordered phase (here h = 1) : the RG flows of the typical 
value (In J^ HP ) of Eq. [28] and of the log of the width ln( Ai n j L ) of Eq. [27] display slope 1 and slope (—1/2) respectively. 



where d s = d — 1 = 1 represents the interface dimension, and where a represents the surface tension. We are able to 
measure the asymptotic behavior of Eq. [28] only sufficiently far h < 4.25 from the critical point h c ~ 5.15, so that we 
cannot measure the critical behavior of the surface tension a(h). 

As shown on Fig. f?l the width Ai n j L of the distribution of the logarithm of the external renormalized coupling 
decays as 

A lnjL cx L-°- 5 (29) 

L— >-\-oo 



V. ANALYSIS OF THE CRITICAL REGION 
A. RG flow of the renormalized transverse field at criticality 




FIG. 5: (a) Statistics of the renormalized transverse field fei at criticality h c ~ 5.15 : (— ln/i^) and Ainh^ as a function 
of L (Inset: log-log plot of the same data) (b) Test of the finite-size scaling form of Eq. [32] : the data collapse obtained with 
ip ~ 0.5 and vfs — 1.3 is satisfactory 



At the Infinite Disorder critical point, the renormalized transverse field is expected to display an activated 
scaling in L with some exponent -0 



ln/i£ oc —L^v c 



(30) 
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where v c is an O(l) random variable. We show on Fig. [5] (a) our data concerning the RG flows of the typical value 
(ln/ii,) and the width A\ n h L °f the distribution of the logarithm of the renormalized transverse field at criticality 
h = h c — 5.15. Our data are consistent with the same scaling of both, with an exponent in the region 

0.4 <ip< 0.5 (31) 

in agreement with previous estimates based on standard Strong Disorder RG [H-Qjl or on quantum Monte-Carlo 
"201 21 1 . The numerical estimate of Eq. [3T1 from the scaling at criticality is not precise as a consequence of the 



uncertainty of the exact location of the critical point (Eq. [T3j) . and of the curvature of the data in log- log plots. It 
is thus interesting to discuss the finite-size scaling in the critical region to relate tj; to other critical exponents of the 
ordered and disordered phases measured in previous sections. 

In the critical region around this Infinite Disorder fixed point, one expects the following finite-size scaling form for 
the typical values p 

ln^ p = Wh£ = -L 4 'F h (L 1 ^ FS \h - h c \) (32) 



where vps is the correlation length exponent that govern all finite-size effects in the critical region. (The exponent 
vps is expected [H, H| to correspond to the exponent v av of the averaged two-point correlation function vps — v av ). 

The compatibility of the finite-size scaling form of Eq. [32] with the essential singularity of Eq. [22] concerning the 
disordered phase yields the relation 

K = iJjv fs (33) 
whereas the compatibility with the behavior of Eq. [M] concerning the ordered phase yields the relation 

vh = (l - ~) vfs (34) 

Eliminating vps, we may thus obtain a numerical estimate of ip from our previous measures of k ~ 0.65 (Eq I23p and 
Vh ~ 1 (Eq[H 

V = — ^ ^ 0.49 (35) 
and the corresponding finite-size exponent then reads 

vfs = v h + %* 1-32 (36) 
a 

These two values are close to the values ip ~ 0.48 and v ~ 1.25 obtained by standard Strong Disorder RG [Til. Ii~2l. Il4l. 
[l5j . To test Eq. [32] we show on Fig[5](b) the satisfactory data collapse obtained with tp ~ 0.5 and vps — 1.3. 



B. RG flow of the renormalized external coupling Jl at criticality 

We show on Fig. [5] (a) our data concerning the RG flows of the typical value In J l ^ v and of the width Ai n j L of the 
distribution of the logarithm of the renormalized external coupling at criticality h = h c ~ 5.15. Our data are again 
consistent with the same scaling of both, with an exponent ip again in the interval of Eq. 1311 

To estimate a more precise value, it is again interesting to consider the finite-size scaling properties in the critical 
region Q 

In J* OT = hTjZ = -L^'Fj (L l ' VFS \h - h c \\ (37) 

The compatibility the finite-size scaling form with the behavior of Eq. [15] concerning the disordered phase implies the 
following relation between exponents 

v typ = (1 - ip)v FS (38) 

The comparison with the relation of Eq. [33] yields with our previous measures of v tyv ~ 0.64 of Eq. [T7] and of 
k ~ 0.65 (Eq[23]) 



(39) 
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FIG. 6: (a) Statistics of the renormalized external coupling Jl at criticality h c ~ 5.15: (— In J^ yp ) and Ai n j £ as a function 
of L (Inset: log-log plot of the same data) (b)Test of the finite-size scaling form of Eq. [37J : the data collapse obtained with 
ij) ~ 0.5 and vps — 1.3 is satisfactory. 

in agreement with Eq. 1351 The corresponding value 

VFS = Vtyp + K ~ 1.29 (40) 

is close to the value of Eq. 135] and to to value v ~ 1.25 obtained by standard Strong Disorder RG [ll|, [IH, [HI, [l5j]. To 
test Eq. (37J we show on Fig|5](b) the satisfactory data collapse obtained with ip — 0.5 and vps — 1*3. 

As a final remark, we have thus found that the activated exponent ip ~ 0.5 at criticality is greater than the 
fluctuation exponent uj = 1/3 of the disordered phase (Eq. [121), as already found on fractal lattices in d > 1 [l8[, 
whereas in dimension d = 1, the two exponents coincide ip{d = 1) = 1/2 = uj(d = 1). This means that the amplitude 
A(h) of Eq. IT4l should diverge as A(h) oc (h — h c )~^~^ VFS , but we are not able to measure this singularity with our 
data. 



VI. CONCLUSION 



To avoid the complicated topology of surviving clusters induced by standard Strong Disorder RG in dimension 
d > 1, we have introduced a modified procedure called 'Boundary Strong Disorder RG' for the Random Transverse 
Field Ising model in d — 2 and we have presented the corresponding numerical results. We have found that the 
location of the critical point, the activated exponent -0 ~ 0.5 of the Infinite Disorder scaling, and the finite-size 
correlation exponent vpg ~ 1.3 are compatible with the values obtained previously by standard Strong Disorder RG 
We thus believe that our modified simplified procedure captures correctly the critical properties. In addition, 
we have analyzed in detail the RG flows within the two phases. Within the disordered phase, we have measured the 
typical correlation exponent vt yv ~ 0.64, which is very close to the correlation exponent Vp Ure {d = 2) ~ 0.63 of the 
pure two-dimensional quantum Ising Model; the fluctuation exponent uj ~ 0.35 which is compatible with the Directed 
Polymer exponent ujdp(1 + 1) = | in (1 + 1) dimensions, in agreement with the arguments of Ref [22j |; the essential 
singularity exponent k ~ 0.65. Within the ordered phase, we have measured the typical exponent ~ 1, which is 
close to the value Vh{d = 1) = 1 in d = 1. 

The simple values found here for v tyv and Vh in d = 2, together with the exact solution in d = 1, raise the question 
whether, in higher dimensions d > 3, the typical exponent v typ still coincides with the correlation exponent of the 
pure quantum transition v tyv (d > 3) = v^ ure {d > 3) = 1/2 ? and whether the typical exponent Vh still keeps the 
simple value vu{d > 3) ~ 1 ? 

More generally, we hope that the idea to change the order of decimations with respect to the standard Strong 
Disorder RG in order simplify the renormalized spatial structure will be useful for all types of models controlled by 
Infinite Disorder scaling. 
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Appendix A: Reminder on Strong Disorder RG rules on arbitrary lattices 

For the Random Transverse Field Ising Model of Eq. [3l we recall that the standard Strong Disorder Renormalization 
are formulated on arbitrary lattices as follows Q : 

(0) Find the maximal value among all the transverse fields hi and all the ferromagnetic couplings Jjk 

ft = max [hi, Jjk] (Al) 

i) If D, — hi, then the site % is decimated and disappears, while all couples (j, k) of neighbors of i are now linked via 
the renormalized ferromagnetic coupling 

J jk = J jk H — (A2) 

ii) If n = Jij, then the site j is merged with the site i. The new renormalized site i has a reduced renormalized 
transverse field 

h ne W = hp_ (A3) 
Jij 

and is connected to other sites via the renormalized couplings 

j-r = ^ + j 3 k (A4) 

(hi) return to (0). 

These standard Strong Disorder RG rules should be compared with the modified procedure called 'Boundary Strong 
Disorder RG' introduced in section [TT1 
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